SUBMERSIONS, HAMILTONIAN SYSTEMS AND 
OPTIMAL SOLUTIONS TO THE ROLLING MANIFOLDS 

PROBLEM 



ERLEND GRONG 



Abstract. Given a submersion tt : Q — )■ A/ with an Ehresmann con- 
nection T-L, we describe how to solve Hamiltonian systems on M by hfting 
our problem to Q and doing our computations here. We also show how 
solutions of lifted Hamiltonian systems on Q can be viewed as particles 
in M under the influence of a generalization of the Lorentz force. We 
investigate what this means in terms of normal and abnormal extremals 
in optimal control problems on respectively M and Q. In particular, we 
are interested in submersions between sub-Riemannian manifolds. We 
apply our theory to a concrete example, namely solving the problem 
of how we can roll one Riemannian manifold M on another Af without 
twisting or slipping along curves of minimal length, connecting two given 
configuations. Here, we both use a submersion to lift our problem to 
a space where computations are simpler to execute, and we show that 
the curves in each of the manifolds can be considered as trajectories 
of objects with a non-constant gauge-charge under the influence of the 
Lorentz force. The Yang-Mills field is in this case given by the difference 
of the Riemannian curvatures of the manifolds M and M. In particular, 
we show that in the case of two surfaces rolling on each other, solutions 
are determined by the equation of a pendulum whose length is inverse 
proportional to the difference in Gaussian curvature of the two surfaces. 



1. Introduction 

Consider a submersion tt : Q ^ M between two Riemannian manifolds 
satisfying the following condition. Let V be the subbundle of TQ formed by 
the kernel of the differential vr* and use Ti for its orthogonal complement. 
Then vr is called a Riemannian submersion if vr* is an isometry on each fiber 
when restricted to Ti. Several papers exist connecting the geometry of Q 
with M, see e.g. [HI El US]- One of the possibilities such a submersion gives 
us, is that we are able to find the geodesies of M, by looking for geodesies 
in Q which are horizontal to Ti at one point (and hence all points) and then 
projecting them to M. This can be a great advantage if it is easier to do 
computations on the space Q than it is on M. A simple example of this, is 
the case the case when Q is a Lie group and M is a symmetric space. 
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A similar idea was used in |12j , for the case when the top space is a sub- 
Riemannian manifold. For any Ehresmann connection 7^ on a submersion 
TT : Q M into a Riemannian manifold M, if we give Q a sub-Riemannian 
structure by lifting the metric of M to Ti, then geodesies in M lift to nor- 
mal geodesies in Q. For the definition of sub-Riemannian manifolds and 
normal geodesies, see Section 13.31 Both Riemannian geodesies and normal 
sub-Riemannian geodesies may be considered as solutions to Hamiltonian 
systems. We will show that given any Hamiltonian function H on M and 
an a choice on Ehresmann connection T-L on the submersion tt : Q ^ M, we 
can define a lifted Hamiltonian function H on Q. We then describe how we 
can obtain the solutions H on M by looking at some of the solutions of 
H on Q. Using this result, we are able to solve any Hamiltonian system on 
M by doing computations on Q. 

We also show a converse result inspired by the following example. Con- 
sider a principle G-bundle n : Q ^ M, where the base manifold M is a 
Riemannian manifold. Again, let H be an Ehresmann connection on tt, 
invariant under the action of G, and lift the metric of M to H, giving Q 
the structure of a sub-Riemannian manifold. In this case, projections the 
normal geodesies of Q are the trajectories of gauge-charged particles under 
the influence of the Lorentz force [13]. Here, the gauge is an element in 
the dual Lie algebra of G and the Yang Mills field is represented by the 
curvature of H. See Section [3l4l for details. Generelizing this idea, we show 
hos solutions of lifted Hamiltonian systems on Q can be viewed from the 
base space M. These projections are curves following a combination of the 
Hamiltonian vector field of the original Hamiltonian, along with a "force" 
depending on the curvature of Ti. The concept of gauge is replaced by an 
introduction of parallel transport for forms which vanish on H. By doing 
these modification, we are able to describe what way solutions are influenced 
by the original Hamiltonian and the connection T-L, respectively. 

The paper is structured as follows. In Section [2] we describe how to lift a 
Hamiltonian system from the base space to the top space of a submersion. 
Then we present our main theorem, connecting solutions of the two Hamil- 
tonian system. In Section[3]we apply this result to optimal control problems, 
in particular sub-Riemannian manifolds. Given a submersion vr : Q — )■ M, 
we both describe how to solve an optimal control problem given on M by do- 
ing computations on Q, and how solutions of lifted optimal control problems 
on Q look when viewed from M. 

In Section H] we give a specific example to show the effectiveness of our 
approach. We look at the problem of finding the optimal curves of the kine- 
matic system of two Riemannian manifolds M and M rolling on each other 
under the constraint of high friction, preventing slipping and twisting. Here, 
results has, up until now, only been presented in some specific cases [TOlfTT]. 
and even the case of surfaces with non constant curvature rolling on each 
other was unknown. We present the equations solving this optimal control 
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problem using exactly the techniques described above. First, we lift the 
problem from our original configuration space to the product of the orthog- 
onal frame bundles F{M) x F{M), where it is simpler to do computations. 
Then, we describe the solutions by looking at the curve projected to one 
of the manifolds. We find that these projected curves behave as objects 
with an 5o(n)*-gauge under the the influence of the Lorentz force, where 
this gauge changes according to a given vector field. The Yang-Mills field 
is in this case given by the differences in the Riemannian curvatures. In 
particular, we show in the case of two surfaces rolling on each other without 
twisting or slipping, the solution of this problem are in many cases found by 
solving the equation of a pendulum, whose length is inverse proportional to 
the differences in Gaussian curvature of M and M along the curve. 

Some of the proofs, including the proof of the main theorem, are left to 
Section [H 

2. Lifting Hamiltonian systems 

2.1. Notation and converntions. For any vector bundle : E — ?> M, 
we will use T{E) is the space of smooth sections of E. Rather than writing 
X{m) for the value of the section X:M— )-i?atmEM, we will use the 
notation X\m- The zero section of the vector bundle will be denoted by 0^. 
If ii^ is a subbundle of the tangent bundle, then a curve 7 in M is called 
£^-horizontal if it is absolutely continuous and has its tangent vector in E for 
almost every t. For elements 61,62 G E^ in the fiber over the same element 
m G M, define the vector vlg^ 62 G Tf.^E by 

d 



vlei 62 



dt 



(ei + 162). 
=0 



This is called the vertical lift of 62 to 61. Similarly, for any X G r(£'), we 
can define the vertical lift vlX G T{TE) of X by vlX|e = vie X\YiE(^^y 

Rather than introducing a special notation for the interior product, we 
will for example denote the interior product of a two-form by a vector field 
X by •). Its pullback by a map / will be denoted r]{X, •). 

2.2. Ehresmann connections, curvature and parallel transport. Let 

vr : Q ^ M be a submersion between two connected manifolds Q and M. By 
a submersion, we mean a surjective map such that the differential map vr* is 
surjective as well. We assume that the manifolds Q and M have respective 
dimensions n + u and n. For each m G M, we write Qm = 7r"^(m) for 
the preimage of m. In what follows, objects related to Q will typically be 
marked by a tilde (~). 

We call the subbundle V = kervr* the vertical space. An Ehresmann 
connection ?^ on vr : Q — > M is a subbundle of TQ satisfying TQ = HQV. 
Since for every q G Q, the map vr* j-^^ : Tig — )• T^(^q^M is a linear isomorphism, 
a choice of Ehresmann connection allows us to define horizontal lifts of 
vectors, vector fields and curves. 
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To be more specific, for every vector v G T^M and q G Qm, we define the 
horizontal hft hqV of v to q as the unique element in Hq which is sent to v 
by TT*. Similarly, any vector field X G r(TM) has a horizontal lift hX given 
by formula 

hX\q := 

Finally, a curve 7 in Q is a horizontal lift of an absolutely continuous curve 7 
in M if 7 is ^{-horizontal and satisfy it{'j) = 7. In other words, 7 is a solution 
of the equation 7 = h^^. Clearly this implies that 7 is uniquely determined 
by 7 up to initial condition. For a general submersion tt : Q — > M and curve 
7 in M, we only know that a horizontal lift exist for short time. 
The curvature operator corresponding to H is defined as 

Ti{X, Y) = pry [pr^ X, pr^ y] for any X,Y e r{TQ) . 

Here pry and pr^ are respective projections to T-L and V, having the other 
bundle as their kernel. This gives us an operator TZ : T(TM) x r(TM) — >■ 
r(V) defined as n{X,Y) = n{hX,hY). It is simple to verify that 7^ is 
(M)-linear in both arguments. 
We also want to introduce an idea of parallel transport of sections in V. 
We do this by defining an operator 

T : r(TM) X r(v) ^ r(v) 

(X,T) ^ TxT:= [hX,T] ' 

This is M-lincar in the second argument and C°°(Af )-linear in the first ar- 
gument. Hence only depends on -'^^(f;) in the first argument. 

Lemma 1. Let 7 : [0, ti] M he an absolutely continuous curve with 
7(0) = mo- Let go be a given point in Qmo, o,nd let 7 he the horizontal lift of 
7 with 7(0) = go- Then for any Tq G Vq^^, there is a unique vertical vector 
field T{t) along 7 with T^(j)T(t) = and T(0) = To. 

Proof. Let {x, U) be a chart on M around niQ, and chose a neighborhood U in 
7r~^{U) around go sufficiently small such that V \ fj is trivial. Let Ti, . . . , T^, 
be sections of V, forming a basis on U. Write 

and 7(t) = X^^^^ Xi{t)dxi. Then T{t) = J2k=i ^«(*) ^kItW ^ solution to 

n u j n u \ 

These equations along with the initial conditions uniquely determines T. □ 

Remark 1 . In contrary to parallel transport with respect to an affine connec- 
tion in a vector bundle, we are only sure that we are able to define parallel 
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transport for short time, since the horizontal hft of 7 may not be defined 
for aU time. 

2.3. Induced submersion of the cotangent bundle and lifted Hamil- 
tonian system. Given a choice of sphtting TQ = 7^ © V of TQ, we will 
also have a corresponding splitting of the cotangent bundle T*Q = Ann(V)© 
Ann(^). By Ann(V) and Ann(^) we mean elements of T*Q which vanish 
on respectively V or 7^. Note that for each q & Q, we have a linear isomor- 
phism T*^^^M — )• Ann(V)g given by p 1— )• vr*(p). Corresponding to our choice 
of subbundle H, we have a unique way to extend the inverse of this map to 
a submersion of the entire cotangent bundle. 
Let us introduce the vector bundle map 

T*Q-^T*M 

YlT*Q jjT*M 

Q 

uniquely defined by the properties kervr^ = Ann(?^) and 7r'^{p){v) = p{hqv) 
for any g € Q, p G T*Q and v € T*^g^M. 

With this formalism, we are now ready to present our main theorem. Let 
H € C°°{T*M) be a Hamiltonian function, and use H for its Hamilton- 
ian vector field. By abusing terminology slightly, we will often refer to an 
integral curve X{t) of H as simply an integral curve of H. 

Theorem 1. Let H G C°°{T*M) and define H G C°°{T*Q) as H = Ho-k'^. 

(a) A curve in X{t) is an integral curve of H if and only if it is, at least 
for short time, the projection of an integral curve X{t) of H which 
is contained in Ann(V) (it is sufficient to require this in only one 
point). 

(b) LetX : [0,ti\ T*Q be an integral curve ofH withj{t) := n'^*'3(A(t)) 
and vo := A(0) pry. Let 7r'^(X) = X with 7 = n^**^(A). 

Then X is a solution to 

(1) A(t) = ^U(t)-vU(j)t;(t)7l(7(t),-) 

where the curve v{t) in Ann{'H) determined by equations 

(2) w^^t)v{t) = 0, v{0) = vo. 

The curve 7 can be obtained again as an H-horizontal lift of ^. 

In the above theorem, v{t)TZ{'y{t), ■) £ T*^^-^M is the element such that 
for any Yq £ T^^i^M it takes the value v{t)TZ{X,Y), where X and Y are 
any vector fields with = 7(i) and Y\^(^)--^ = Yq. Parallel transport in 

(I2j) is defined such that for any X G T{TM), T G r(V) and a G r(Ann(7^)), 
we have {Wx Oi){T) = X{a{T)) — a{Wx T")- 

The proof of this theorem is left to Section 15. 1[ 
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Remark 2. In this presentation, we started by choosing an Ehresmann H 
and used this to construct vr^ and H. Conversely, suppose that we have a 
Hamiltonian H which is equal to H o ■k'^ for some H € C{T*M) and vector 
bundle morphism tt'^ : T* Q ^ T* M over vr which is a left inverse to vr* on 
every fiber. Then we can define Ti as the vectors annihilated by kerrr^, and 
the same result holds. 

2.4. Special case: Hamiltonian of a Riemannian manifold. We now 

consider the following special case. Let (M, g) be a Riemannian manifold 
with metric g. Associated to this metric, let b : TM — >■ T*M be the bundle 
isomorphism v >—?■ g(f , •) with inverse ft- Use V for the Levi-Civita connection 
corresponding to g. Associated to the Riemannian structure, we also have 
a Hamiltonian function H{p) = ^p(jjp). The projection of integral curves of 
H to M gives us the Riemannian geodesies. If we apply Theorem ^h) to 
this Hamiltonian we get the following result. 

Corollary 1. Let ir : Q ^ M be a submersion into a Riemannian man- 
ifold (M, g), with associated Hamiltonian H. Let H C TQ be a chosen 
Ehresmann connection, and define H as in Section \2. 3X 

Let A : [0, ti] — ?> T*Q be an integral curve of H with A(0)pry = vq. Use 7 
for its projection to Q. Then ^ is a horizontal lift of a curve 7 in M , which 
is a solution of 

v^it)i{t) = -^v{t)iz{^,-), 

where W-y v{t) = and v{0) = vq. 

Proof. Pick a local basis ei, . . . , en be a local basis orthonormal vector fields, 
and use these vector field to give the cotangent bundle local coordinates 
Pi = Pi^i)- It is standard result that the Hamiltonian vector field in these 
local coordinates is given by 

n n 

H = ^Piei- ^ pjPkT'ji^dp^, r^^. = g(efe,Ve,ej). 

i=l i,j,k=l 

Let A(t) = 7r^(A(t)). If A(i) = J2i'=iPi{t)'<^^i\'y{t) is a solution to 

m = H\^^t) -^h(t)v{t)n{j{t),-), 

thenbe,(7(t)) =K(t) andp,(t) = - E",fc=iPi WPfc Wrjfc^W-^W ^(^W' e*)- 
The corollary then follows from the computation 



n n 

'^^{t)i{t) =^Pi{t)ei+ PiPj^ijek 

i=l i,j,k=l 
n j n \ n 

= Y.\-Y. P3Pk^]k - <t) ^(7(i), e.) + ^ PiPjT^eu 

i=l \ j,k=l I i,j,k=l 



= -Ht)mit),-). 

□ 
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3. Lifted optimal control problems 

3.1. Optimal control systems. Then are many different definitions and 
generalizations of an optimal control problem. We will use the definition 
found in [H Section 2.1] and [2]. 

A smooth control system consists of a fiber bundle ^ :IA ^ M with fiber 
U , along with a bundle morphism 

U ^TM . 




A measurable essentially bound curve K{t) inU is called an admissible control 
if its projection 'y{t) in M is an L°° curve satisfying 7 = f{K{t)). We want 
to consider the following optimal control problem: For a smooth function 
if : C°°{U) M and two points mQ,mi S M, find the admissible control 
K{t) which satisfies 

7(0) = mo, 7(ti) = mi, with = tW, 

and minimize the functional 

I Lp{K,{t)) dt. 

Jo 

A sufficient condition for an admissible control to be a solution to this opti- 
mal control problem is given by the Pontryagin Maximum Principle (PMP). 
In order to present this result in a simpler way, we will just write the formu- 
lation assuming that U can be written as a trivial fiber bundle U = M x U , 
which can be considered a local version of the general case. See Remark [3] 
for how this statement can be reformulated for the case when U cannot be 
trivialized. 

Theorem 2. PMP for Optimal Control Problem with fixed time ti 

[II Theorem 12.3]: 

Let K{t) = {j{t),u{t)) be a solution to the above optimal control problem, 
where j{t) is a curve in M and u{t) is a curve in U . For each r (^M,u (z U , 
consider a Hamiltonian function 

(3) Hr,u{p) = p{f{rn,u)) +rip{m,u), p€T^M. 

Then there exists a curve A : [0, ti] — ?> T*M, and a number r < such that 

(i) n^*^(A(t))=7(t). 

,u{t)\\{t) for almost every t, 
(iii) Hj.:^n)(X(t)) = max Hr u(Mt)) for almost every t. 

Moreover, if r = 0, then A never intersects the zero section ofT*M. 

Let us use the name extremals for solutions of PMP. They are called 
normal if r 7^ (it is sufficient to consider r = — 1) and abnormal if r = 0. 
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Remark that abnormal extremals to not depend on the function ip, only the 
control system (Cf)- 

Remark 3. If the fiber bundle cannot be trivialized, PMP can be reformu- 
lated in the following way. For any r G M, define J^r ■ U x mT*M R 
by J^r{n,p) = + rip{K,). which takes the place of the Hamiltonian 

in ([3]). Requirement (ii) is then replaced with the identity 

(pi'2?)k{t),A{t)(^, A) = (Jf r)*|{K(t),A{t))(^)- 

for any vector field X G T{T{U x mT*M)), where pr2 : U xmT*M T*M 
is the projection and <; is the canonical symplectic form on T*M. In the 
requirement (iii) , the maximum needs to hold over all elements in Z^^(t) • 

3.2. Submersions and a liftet optimal control problem. Let the pair 
^ -.U ^ M and f : U ^ TM be a control system and let vr : Q — ;> M be 
a submersion. Then a choice of Ehresmann connection % on vr gives us a 
unique way induce a control system on Q. 

• Define a fiber bundle ^ : vr* — )■ Q as the pull-back of ^ : U ^ M. 
Recall that this is defined as 

TT*U = {iq,K) eQ xU : TT{q)=C{K)}. 

The map ^ is just the projection on the first coordinate. Let vf denote 
the projection on the second coordinate. Then we have the following 
commutative diagram 



• Define a bundle morphism f : tt*1/( ^ TQ by 

/(k) = for any Ken* U . 

In other words, if k represents the pair {q, k), then /(k) = hqf{K,). 

On this control system, we consider the optimal control problem induced 
by the function if = (p on. In this situation, extremals of Q and M, both 
normal and abnormal, are connected through Theorem [TJ It is sufficient to 
show this locally, so we assume that we can identify vr* lA with Q x U. We 
then verify that for any p £ T*Q and u £ U, we have 

Hr,u{p) ■= pU{(Ii u) + r!p{q, u), 

= p{hqf{TT{q),u)) + rip{TT{q),u) 

= ^\p){f{<<l).u)) + Mn{q),u) = Hr,u{^\p)). 

This relation can be a powerful and useful tool in many optimal control 
problems. We will explore the case of sub-Riemannian manifolds in detail. 
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3.3. Submersions and sub-Riemannian manifolds. A sub-Riemannian 
manifold is a triple (M, D, g), where M is a connected manifold, Z) is a sub- 
bundle of TM and g is a metric tensor on D. Use AC^) for the Hilbert 
manifolds of all D-horizontal curves defined on the interval [0, 1] with an 
derivative. Let ACi:)(mo, rrai) be the subset of ACd consisting of curves7 
satisfying 7(0) = rriQ and 7(1) = mi. The distance function dec in this 
geometry is called the Carnot-Caratheodory distance and is defined as 

dcc{'mo,mi) = inf i^j g(7, 7)^/^ : 7 G ACD(mo, mi)| . 

A sufficient condition for this distance to be finite, i.e. that any pair of points 
can be connected by a D-horizontal curve, is that D is bracket generating. 
This means that D along with the iterated brackets of its sections span the 
entire tangent bundle. To look for minimizers with respect to this distance 
function, is equivalent to look for curves in ACdCw-O; w-i) which minimize 
the energy functional 

1 

E{l) = ^J^ dt. 

This can be viewed as the control system with data (n^,inc), such that 
n-^ : D — )• M is the projection and inc : D — ?> TM is just the inclusion, 
where we try to minimize £'(7). The only difference is that we allow the 
derivatives to be in rather than L°°J3 

Normal extremals can be described in the following way. The metric ten- 
sor g defines a bundle map jjg : T*M D hy the identity g(tigP, v) = 
p{v) for any v € D. Normal extremals are then the solution of the Hamil- 
tonian system with Hamiltonian Hsr{p) = ^p{^gp). Projections to M of 
normal extremals are always smooth and are length minimizers locally. We 
therefore call such curves in M normal geodesies. Remark that if A is a 
normal extremal, and its projection to M is the normal geodesic 7, then we 
have the relation 

(4) at) = km 

Abnormal extremals A also have an alternate description. Let <j be the 
canonical symplectic form on T*M and let A be a curve in T*M. Then A 
is an abnormal extremal if and only if it is an absolutely continuous curve 
in Ann(D) with an derivative, never meeting the zero section, such that 
?(A(t), •)|Ann(D) = 0. A curve fulfilling the latter requirement, is often re- 
ferred to as a characteristic of D. We will use the term an abnormal curve 
7(t) if the curve is the projection to M of an abnormal extremal X{t). Ab- 
normal curves are singular points of the map from the Hilbert manifold 
AC£)(mo) of curves in ACd starting at mo to the manifold M given by 

end : ACD(mo) ^ M, end(7) = 7(1). 

^Some authors prefer to work with sub-Riemannian geometry, requiring that horizontal 
curves should have L°° derivatives. 
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For more details on sub-Riemannian manifolds, see 

Let TT : Q — )• M be a submersion, with vertical space V and a chosen 
Ehresmann connection Ti. We lift the sub-Riemannian structure in the 
same way that we do with a more general optimal control problems. The 
lifted structure can be described as the sub-Riemannian manifold {Q,D,g), 
where 

D = {hgv : V e D^,Q eQm,me M} = {TT^)-\D)nn, 

and g = vr* g . We will say that this sub-Riemannian structure on Q is lifted 
from M. We present the following corollaries of Theorem [TJ 

Corollary 2. A curve X{t) in T*M is a normal extremal if and only if it is 
a projection of a normal extremal in T*Q contained in Ann(V). Conversely, 
projections of normal extremals in T*Q satisfy equation ([T]) with H = HgR. 

Remark 4. In the special case when D = TM, making the base space a 
Riemannian manifold, projections of normal geodesies are given by Corol- 
larylll The top space is then a sub-Riemannian manifold (Q, H, vr* g). Also, 
^-horizontal lifts of Riemannian geodesies are normal geodesies. The latter 
fact was first observed in [12, Th. 6.2, Cor 6.5] 

Theorem 3. A curve 'y in M is abnormal if and only if any horizontal 
lift of the curve is abnormal. Conversely, an L'^ -curve X{t) in T*Q with 
TT^m) = X{t), U^'Q(X{t))=m, v{tl=m piv andW^'^^iXit)) = 7(t), 
is an abnormal extremal if and only if j is D -horizontal, 

<^{K-)u^^0) = -v{t)nj,nr^'-), 

and A and v does not vanish simultaniously. 

The proof of this theorem uses elements of the proof of Theorem [H and 
is therefore left to Section [ 



3.4. Special case: Charged particles in a magnetic field or a Yang- 
Mills field. Consider the special case when the submersion tt : Q — )■ M is 
a principal G-bundle, where G acts by a right action. Corresponding to this 
action, we have vector fields cr(A) defined such that for any A in the Lie 
algebra g of G, 



(5) = I 



q ■ exp(j{tA). 

=0 



For every q, the map i— V |g defined by A i— > cr(^)|g is a linear isomorphism. 

Let Ti be an Ehresmann connection which is invariant under the group 
action, that is. Hp -a = Tip. a for any a £ G. Such an Ehresmann connection 
is called principal, and corresponding to it is a principal connection form uj, 
which is a g- valued one-form, defined by having T-L as its kernel and satisfying 
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uj{(T{A)\q) = A. It is simple to verify that if T is a vertical vector field along 
an "H-horizontal curve 7 in Q, and tt(^) = 7, then 

T = 0, if and only if T(t) = a{A)\^(^f) for some A € . 

Introduce the curvature two-form 0, by 0,(X,Y) = du){X,Y) + a;([X,y]). 
Then we can rewrite ([1]) as 

X = H\x-ylxLn{j,-), 

for some constant L € 0*. If M is a Riemannian manifold, and H is the 
corresponding Hamiltonian, then we can write Corollary [1] as 

V^7 = -tJLfi(7,-). 

These are the equations for a particle with gauge L under the influence of 
the Yang-Mills field O. If G = U{1) or M, then L represents the charge and 
n is a magnetic field. For more details, see [13] and ^14^ Chapter 12]. 

4. Optimal control of rolling manifolds 

4.1. Definition of manifolds rolling without twisting or slipping. In 

this section, unless otherwise stated, M and M will denote connected ori- 
ented Riemannian manifolds, both of dimension n. The Riemannian metrics 
of the respective manifolds will be denoted by g and g. To avoid trivial con- 
siderations, we will always assume that n >2. 

We consider the kinematic system of M rolling on M. The configuration 
space for this motion can be described as follows. Define the fiber bundle Q 
over M X M as 

(6) Q = {sO(r^M, T^M) : m G M, m G it?} . 

Here, SO(V^, V) denotes the space of all linear, orientation preserving isome- 
tries between two oriented inner product spaces V and V. Clearly each fiber 
is diffeomorphic to SO(n). An element q : TmM — > T^M in Q represents 
a configuration where M lie tangent to M at the points m and m. Two 
vectors, one in each tangent space, lie adjacent if one is mapped to the other 
by q. A rolling is then a curve in this configuration space. 

We will assume that we have high friction, giving us the constraints that 
we cannot slip or twist. By slipping, we mean moving one of the manifolds 
while remaining connected in the same point on the other. By twisting, we 
mean spinning in place, that is, we move in the same fiber over M x M, 
changing configuration, but not connecting points. Mathematically, we can 
write these constraints as follows 

Let TT : Q — )■ M and vf : Q — )■ M be the respective natural projections. 

Definition 1. An absolutely continuous curve q{t) in Q, with ^{t) = Ti{q{t)) 
andj{t) = 7f(g(f)), is a rolling without slipping or twisting, if it satisfies the 
following conditions for almost every t, 
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• (No slipping condition:) q{t)j{t) = 7(t), 

• (No twisting condition:) Any vector field X{t) along j(t) is parallel 
if and only if q{t)X{t) is parallel along 7(t). 

In what follows, we will typically drop the phrase "without twisting or 
slipping" and just refer to a rolling, with the constraints being implicit. For 
more details on these types of systems, see e,g. [151 El IS]. 

Rolling without slipping or twisting can be described as horizontal curve 
relative to a subbundle D of TQ of rank n. We describe this subbundle lo- 
cally. On a sufficiently small neighborhood U on M, choose an orthonormal 
basis ei, . . . , e„ of vector fields on U . Define U and ei, . . . , similarly on 
M. We can then trivialize the Q over [/ x [/ by 



Q\u^u ^ C/x;7xS0(n) 

qeSO{TraM,Ti^M) ^ {m,m,{qij)) ' 

Relative to this trivializations, define vector fields 



Qij = siei\fn,q{t)ej\ 



s=l 

Then D on Q]^^^ is spanned by vector fields 

(7) ej = ej + qej+ ^ (g(eQ, Ve^e/j) - g(gea, Vge^ge^)) 

l<a</3<n 

Here, qcj denotes the vector field q i-> (J'e|^(q). 

Notice that D is an Ehresmann connection for both tt : Q M and 
vf : Q — )■ M. We can lift the Riemannian metric g on M to a metric h on 
D, by 

h{vi,V2) = g{Tr^Vi,-K^Vi), Vi,V2eD. 

In this metric, the vector fields in ([7]) form a local orthonormal basis. No- 
tice that from the definition of D, the metric h also satisfies h.(vi,V2) = 
g(7f*wi, 7r*f2)- Hence, the sub-Riemannian structure {D,h) can be consid- 
ered lifted either from M or from M. 

To find a sub-Riemannian length minimizer on Q relative to h with initial 
point qo and final point qi, is to find the rolling q{t) from qo to qi such that 
7(t) = 7r(g(t)) (or equivalently 7(t) = Tf(q{t))) have minimal length. This 
problem can be quite complicated to attack, and few results exists. The first 
results was given for a sphere rolling on a plane [10]. It was shown that in 
this case, projection of geodesies to either the sphere or the plane, are curves 
with constant speed whose direction is a solution of the pendulum equation. 
Later results for more general manifolds of constant curvature are found in 
[llj . In addition, it was shown in [7] that if the Riemannian curvatures R 
and R of respectively M and M satisfy the condition that the map 

f 1 A U2 R{TT^,Vl,7r^V2,TT* VT* •) — i?(7f *f 1 , 7?=,, ^2 , vf* -,7?* •), fl, ^2 G D, 
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from /\^ D to /\^ D* is invertible, then we both have that D is bracket 
generating and there are no abnormal minimizers which are not also normal 
minimizers. Hence, in this case, we can focus on normal geodesies. 

We would like to describe the geodesies in the general case. Unfortu- 
nately, computations in Q defined as in ([6]) has proved to be quite difficult. 
We therefore use Theorem [1] (a) to lift our problem to a space where com- 
putations can be done more easily. We start with some definitions. 

4.2. Frame bundles. We make the convention that whenever we mention 
M", it will always come equipped with the standard orientation and the 
Euclidean inner product. We will write GL(y, V) for the space of all linear 
isomorphisms of vector spaces V to V. A frame at the point m € M, is a 
map / G GL(M", T^M). This can be considered as a basis {/i, . . . ,/„} of 
TjnM, by identifying / with the vectors 

/, :=/(0,...,0,l,0,...0). 

^ V ' 

1 in the j-th place 

Conversely, any choice of basis determines a map /. If we denote by 
J-^(M) = GL(M",r„jM), we can define a frame bundle r : -F(M) M 
as the principal GL(n)-bundle with fibers Trn{M). 

For a given affine connection V on M, we define a corresponding Ehres- 
mann f on r consisting of all tangent vectors of curves f{t) such that each 
fi{t), . . . , fn{t) is parallel vector along Tf{t). This is sub-bundle is invariant 
under the right action of GL(n). 

For an oriented Riemannian manifold, define the oriented orthonormal 
frame bundle r : F{M) — )■ M as the principle SO(n)-bundle, such that the 
fiber over m is Fm{M) := SO(M",rmM), that can be identified with the 
space of all positively orthonormal frames of T^M. 

Let us consider the case where V as the Levi-Civita connection on M. 
Since parallel translation along a curve preserves orthogonality, we can in- 
duce a principal Ehresmann connection £" on r : F{M) — > M instead of 
T{M). This Ehresmann connection again corresponds to an so(n)-valued 
principal connection one-form uj = {ujap) on F{M) defined as in Section [3.41 

Let b : TM — >■ T*M be the vector bundle isomorphism v i— >■ g(v, ■) and 
write jj = b~^. Define an valued one-form 9 = {6i) on F{M), by 

e^\f=r*\>fr 

The differential of 9 and w are connected through the Cartan equations, 

n 

d9i = — ^ Wjfc A 9k- 

k=l 

n 

diOal3 = — ^ ^ak A Wfc^ -|- Qal3- 
k=l 
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where O = (^a/s) is the so(n)-valued curvature two-form. In this case, this 
is connected to the Riemannian curvature by 

n{f){X,Y) = {R{X,YJp,U)U^p. 

R{X, Y, Z, W) = g((VxVy - Vy Vx - y[x,Y])Z, W). 
We can rewrite the Cartan equations in terms of vector fields rather than 
forms. For any A = (Aaf^) G so(n), let cr{A) be defined as in ([5]). These 
satisfy 9i{a{A)) = and a;Q^((T(A)) = A^js. We also have the vector fields 
Xi with values in £, defined by Xi\f = hffi, where the horizontal lift is 
with respect to £. These satisfy 9i{Xj) = 6ij and tOapiXj) = 0. Using this 
information along with the Cartan equations, we obtain 

[XuXj] = -a{n{Xi,X,)). 

n 

\X,,a{A)\ = -Y,A,,X,. 

\a[A\a{B)\=a{\A,B\). 

4.3. Lifting of the rolling problem. To simplify our computations, we 
want to lift our problem to F{M) x F{M). Define a principal SO(n)-bundle 
by 

p : F{M) X F{M) Q 

(/,/) ^ fof-' ■ 

We will use the following result. 

Lemma 2. [3 Cor 1] Let q{t) he a curve in Q which projects to the curve 
(7(^)1 7(^)) in M X M. Then q{t) is a rolling without twisting or slipping if 
and only if it is a projection of a curve (/(t), f{_t)) in Q such that 

(a) fi{t), . . . , /n(t) are parallel along "f{t), 

(b) fi{t), . . . , fn{t) are parallel along 7(t), 

(c) f{t)-'m)) = f{t)-\m), 

Such curves {f{t), f{t)) may be considered as horizontal lifts of the curve 
q{t) with respect to an Ehresmann connection Ti on p. We will describe this 
connection. 

Let 9, u}, 0,, Xj and cr{A) be defined on F{M) as in Section [4.21 and define 
9,Q,i},Xj and ^{A) similarly on F{M). Then 

V = ker p^ = {a{A) + a[A) : A G so(n)} . 
Pick an Ehresmann connection 

v. = span{Xj,Xj,a{A) - cr{A) : j = 1,. . . ,n,A £ so(n)}. 

Use horizontal lifts with respect to H of elements in D to define a distribution 
D. Then D-horizontal curves in F(M) x F{M) are exactly the curves in 
Lemma [2j A basis for D is given by Xj + Xj . Lift the metric h to a metric 
h on D. Then Xj + Xj form an orthonormal basis. By Corollary [21 we 



SUBMERSIONS, HAMILTONIANS AND ROLLING MANIFOLDS 



15 



can now find optimal solutions of the rolling problem by looking for normal 
geodesies of the sub-Riemannian manifold {F{M) x F[M),D,h) which are 
projections of integral curves in Ann(V) of the sub-Riemannian Hamiltonian 
with respect to h. 

4.4. Normal geodesies and optimal solutions. We will use Theorem 
[D^b) to describe the normal geodesies. In order to present our results, let us 
introduce the so(n)-valued two- form on F{M) defined by 

This notation is inspired by the fact that this form coincides with the cur- 
vature form when M = 5". We will also identify so(n) with so(n)* through 
the inner product 

(8) {A,B) = tiiA^B). 

With slight abuse of notation, we also will use b for this identification b : 
so(n) so(n)*, with inverse [j. 

Theorem 4. Let q{t) he a horizontal curve in (Q, D, h) with a projection 
7(t) in M . Then q{t) is a normal geodesic if and only if there is vector fields 
V{t) along j{t) and a curve L{t) in 5o{n)* satisfying 

V^7 = -tt^ (0(7,.)-%7,'?-)) 

«L = -^o^"(7,y), v^v = -m{qi,q-)- 

Remark 5. At first glance, it may seem striking that the curvature tensor 
of appears in the geodesic equations in Theorem SI However, its appear- 
ance is just the result of our choice of inner product in ([8]) to give so(n)* 
convenient coordinates. 

There is also an apparent asymmetry between M and M, since V^l^ 
only depends one the curvature of M. This is of course an illusion. Define 
V{t) = V{t) + 7. Then it is still true that = -\ Vt^"{^, V), but now we 
have V^V = -tti^^(7, •)• 

Proof. Let (/(t), f{t)) be some D-horizontal lift of ^{t) and let q{t) = f{t) o 
f{t)~^. We will use Corollary [1] to find the geodesies, recalling that D is an 
Ehresmann connection on F{M) x F{M) — ?> M. 

We start by finding the curvature IZ oi D. Since elements Xi + Xi form a 
basis, the curvature is TZ{Xi + Xi, Xj + Xj) = a{Q{Xi, Xj)) + a{Q{Xi, Xj)). 

Next, we determine the equations for parallel transport. Write 'y{t) = 
J27=i iii't)^i\'y(t) for some local orthonormal basis ei, . . . , e^j. Introduce local 
coordinates fj = X]i=i fij^i- Using these coordinates 

n 

a{A)= f^sAs,^f^^. 
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Since the D-horizontal lift of Cj is hei = Yl^=i fisi^s + Xg), we have 

We,Xj = -d{Q{hei,Xj)). 

n 

Te, a{A) = - Te, a{A) = ^ /ir-^r^s- 

r,s=l 

If v{t) is a curve in Ann(V)nAnn(L'), then it is on the form Yl^=i 

Oj) + ^ Yl'ij=i —^ij) where j|L = (Lij) is a curve in so(n). In order 

for V to be parallel, it must satisfy 

vj = -dta{Xj) = -a{W^Xj) = -m{h^,Xj) 

n 

{'AL,A) = dta{a{A)) = a(T^cj(A)) = - ^ jifivAsrVs 

i,r,s=l 

Y,{'^if^rVs-'^if^sVr)] , a) = ,V) , A) 




i=l 



where V is the vector field V{t) = Vi(t)fi(t). 

We insert these results into formula Corollary [U and we find 

V^7 = -tt^(i)^(7,-) 



(n 1 " \ 

vjim - + 2 E - •)) + ^mi, •))) 



(0(7,. 



where t|L = -i(17^'(7, F), ^) and V^V = -^Q{q^,q-). □ 

4.5. The two dimensional case. Let us consider the case when the Rie- 
mannian manifolds M and M are two-dimensional. Let fi(t) and f2(t) be 
a basis of orthogonal parallel vector fields along a curve 7(t) in M. Write 

(9) 7 = a(cos6l/i +sin6'/2). 

Use respectively and 0, for the Gaussian curvatures of and M. We will 
make the addition assumption that for any pair of points (m, m) £ M x M 
we have that — ^\fh never vanishes. The rolling distribution is bracket 
generating if and only if this condition holds [HE]. Under this assumption, 
we have the following result. 

Proposition 1. Let q{t)be an absolutely continuous curve in Q, with pro- 
jection 7(t). Write VL{t) = ^\^{t) and Vl{t) = and introduce the 
notation p{t) = ^^^^y^"^^- Represent 7 as in Q. Then q{t) is a normal 
geodesic if and only if a is constant, and 9 is a solution to the equation 
■■ ■ A . 
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where F{t) = p{s)^ sin{e{t) - e{s)){^}{s)n{s) -n{s)n{s))ds and A and (po 
are arbitraty constants. 

Remark 6. If il.{t)n{t) = il.{t)n{t), i.e. if one of the Gaussian curvatures is 
a constant times the other, the geodesic equations reduces to the equation 
a pendulum. More specifically, when F{t) = 0, (jlOp is the equation for the 
pendulum whose length varies by p{t) under the influence of the force of 
gravity with magnitude A and direction 6q. See Section [5^3] for details. 

In particular, this happens if both of the manifolds have constant Gauss- 
ian curvature or if one of the manifolds is flat. 

Proof. It is simple to verify from Theorem [4] that I7I is always a first inte- 
gral. In the notation of the same theorem, let V = bi{cos9fi + sin 0/2) + 
b2{— sin 0/1 -|- cos6f2). Also write t|L and the curvature forms as 

2[-L J ' 2[-n J ' 2[-n 

We can rewrite equations in Theorem [J] in these variables. 

9 = -^L{n - h), 

L = ^ab2, bi = 9b2, 62 = -Obi - ^aLO. 

This can be reduced to three equations by the identity L = = —2p6, 

which give us the formulation 

(11) h = eb2 b2 = -e ( bi + a^L.] . 

\ Q-nJ 

-2pe - 2pd = \0Lb2. 

By using a change of variable and a standard variation of parameters argu- 
ment, we can solve (1111). We the find that 



b2{t) = Asm{e{t) - 0o) + a / sm(0(t) - e{s)) ] ) ' ds, 

Jo vHs) - "(•sj) 

which completes the proof. □ 

4.6. Rolling on M". Assume that M = M" with the Euclidean metric 
and standard orientation. Then we can identify Q with xF{M), by 
associating q : T^M — )• TyW^ with the pair {y,q~^), where we have used 
the canonical identification of TyM"' with M". We then have the following 
reformulation of Theorem SI 

Corollary 3. Let q{t) be a rolling of M on M". 
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(a) Let 7(t) he its projection to M . Then q(t) is a normal geodesic if 
and only if 

where V-^/V = and V^VF = 7, and Lq is a constant element in 
so(n)*. 

(b) Identify q{t) with the curve {y{t),f{t)) in xF{M). Then this is 
a normal geodesic if and only if 

y = -^Lon{f{t)y{t), ■) + mf{t)v, W - y(o)), f{t)y{t), 

for some v G M" and Lq G so(n)*. 

Proof. Here we have used that V becomes parallel when i7 = 0. Also, since 
V, now is parallel, if W is any other vector field along 7, then dtO.^"^ {W, V) = 
n^'\S/yW, V). The result then follows □ 

Let us do a physical interpretation of these equations in the two dimen- 
sional case. The curvature form on M can be considered as an M-valued 
two form, and we write this as = Bvol, where vol is the Riemannian 
volume form induce by the metric and orientation on M. represents a 
magnetic field with magnitude B, always normal to the manifold. 

Let 7(t) be the trajectory of a charged object moving in M. The initial 
charge is Lq, however this charge changes according to a parallel vector field 
V. Let X be a vector field on M extending V locally. Then the charge at 
time t is given by Lq + J^|^^ ^ -kbX, where * is the Hodge star operator. In 

other words, according to the positive orientation, movement in a direction 
rotated 90° from the direction of X will increase the charge, while movement 
in the opposite directions decreases the charge. 



5. Proofs 

5.1. Proof of Theorem [H Before we are ready to present the proof in 
Section 15.1.31 we need to introduce some concepts, such as lifted vertical 
and horizontal bundles and the vr^-symplectic lift of vectors. These concepts 
and some necessary lemmas are described in Sections 15.1.11 and 15.1.21 with 
a representation in local coordinates included in Section I5.1.4[ 

5.1.1. Lifted horizontal and vertical bundle. Let vr^ : T*Q — t- T*M be the 
submersion defined in Section 12.31 relative to some Ehresmann connection 
?^ on TT : Q — 7> M. Prom the first mentioned submersion, we have vertical 
space = ker vr^*. We will show that this subbundle of TT*Q is symplectic 
with respect to the canonical symplectic form on T*Q and we can use this 
to define an Ehresmann connection Ti^ on vr^. 

Elements in the cotangent bundles of Q and M are denoted by respectively 
p and p. We will use for the Louiville one-form on T*M defined by 
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i){w) = p{Il^* ^w) , w € TpT*M and ? = —d-d for the canonical symplectic 
form. Use / forjhe bundle morphism I : TT*M T*T*M, I{w) = (^{w, ■). 
Define i),'^ and I similarly on Q. Relative these structures, we will introduce 
a notation for what we will call the -symplectic lift. For any w G TpT*M 
and p G (7r^)~^(p), we will use IpW for the element 

ipW = I^^TT^*\pI{w). 

Here, by abusing notation, we have used 7r^*|p for the map T*TM — ^ 
T~TQ,ao I-)- ao(7r^*|p •)• other words, IpW is the unique element in 
TpT*Q satisfying 

^{ipW, u) = <j(u;, 7r^*S), for any u G T-TQ. 
We use this notation in the following lemma. 

Lemma 3. LefH? he the symplectic complement ofV^, i.e. = Ann(V^). 

(a) For anyp G T*Q, the fiber ofH^ is given byH~ = {£pW : w G T^2^T*M}. 

(b) For any p G T*Q and w G T^2(^pjT*M, we have 

(c) andy.^ are symplectic, i.e. = TT*Q. Furthermore, 

nT*Q^V2 = V and U^''^,H'' = n. 

Proof. (a) This is an easy consequence of the fact that for any p eT*Q, 
the map T*2^p^T*M — >■ Ann(V^)p, ao i-^ 7r^*Q;o is a linear isomor- 
phism. 

(b) For any two elements p, G T*^M, it is simple to verify that 

?(vlpao,w) = ao(n"^ ^\w), for any w G TpT*M. 

Observe that for any one-form a on T*M, by definition vr^* vl7r*a = 
via. Hence, for any vector field W on T*M, we have 

a([F*^^W) = -^(W^,vla) = -^{eW,Y\7T*a) = a{Tr^U^'^JW). 

Here, the vector field iW on T*Q is defined in the obvious way. Since 

a and W where arbitrary, the result follows. 

(c) Clearly, n'^''^^V^ = n^*'5_^ ker tt^^ C kerTr* = V. This map has to 
be surjective since it is a right inverse to the injective map 0^*'^* : 
TQ TT*Q which maps V into V^. 

To prove the analogous result for T-L^, observe that for any p G T*Q 

and 5o G Ann{'H)q, we have vlp oq G V^. By definition i9(vlp 5o, w) = 
ao(n'^*'5*w) = for any w G Ti^. Hence, n^*'?*^^ ^ ^ ^^^^^ ^j^^g 
map is surjective by (b). 

Finally, assume that w G {n^nV^)\p. From our previous re- 
sults, this intersection is contained in the kernel of n^**^^,. Hence, 
w = vlp 5o for some So GT*Q which is in the same fiber as p. How- 
ever, since ?'(vlp ckq, •) must annihilate both and and since 
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Ti^ U is mapped surjectively onto TM, the identity ?(vlpao,n) = 
ao(n^ ^ *u) implies that qq = 0. This shows that and are 
transversal. 

□ 

5.1.2. Horizontal lifts with respect to 'H'^. A consequence of Lemma [3] (c), 
is that T-L^ is an Ehresmann connection on vr^ : T*Q T*M. Hence, we 
have two different ways of lifting elements in TT*M to T-L^, namely the tt"^- 
symplectic lift and the horizontal lift. We will denote the latter by h-w. 
The following lemma describes the relation of these two notions. 

Lemma 4. Ifp(z T*Q, the following relations hold. 



n 



T*M 



(b) 



(12) 



= epw + Y\ppn{hqlf*^\w, ■). 



Here, p TZ{hqU. 



T'M 



E T*Q is the covectorv i-^ pTZ{hqIl 



T'M 



W, V 



As a consequence 

iT^<t£pW = w — vlpp TZ{Il^ *w, •) P = 7r^(p). 
Remark 7. Written as a commutative diagram, let p in T*Q satisfy 




T*M 



Lemma [H then states that for any that for any w € TpT*M with v 
li^ *w, we have the relations 



w 



and 



hqV 



hqV 



w — vlpp TZ{v, ■ 

YlT*M 



Observe that vlpP TZ{v, •) only depends on the projection of p to Ann('H). 
Hence, it is simple to see that h~w = £pW whenever p is inn Ann(V). 

Proof of Lemma^ (a) This is immediate from the definition of horizon- 
tal lift and Lemma [3] (b) and (c). 
(b) Let us split the Liouville form on T*Q into two parts d = "i?^ + 'd^ 
where 

^'^[w) = p(pr^ lf*^^w), i}^iw) = p(prv U^'^^^w), w G TpT*Q. 



SUBMERSIONS, HAMILTONIANS AND ROLLING MANIFOLDS 21 

We observe that for any w £ TpT*Q, 

= 7r2(p)(n^*A^7r2,S?)=7r2*(^)(^}) 

which in turn imply ^= vr^*? — di?^. 

Since Ti^ is symplectic, h~w is completely determined by the val- 

— — 2 

ues of I{h'~w) on vectors h'-u, where u G T^2^j^TM. Define TZ as 

~ '^2 

the curvature of . From (a), we know that iC'*^ ^TZ {h-w, h~u) = 

TZ{hgU^*^,w,hgU^'^,u). We calculate 
^{h~w, h~u) = q{w, u) — {h~w, h~u) 

= c;{w,u)+^{n^{hlw,hlu)) 

= <;{w,u) +pn{hqU'^*^^,w,hgU'^'^^,u) 

= ?{ipW, h~u) + ^{vlppTZ{hqlf ^.w, •), h~u) 

which shows (jl2p 

□ 

We will state one of the observations made in the proof of Lemma U] as a 
separate result, since we will need it later. 

Lemma 5. If 'd^\p = p(pry n^**^* •) is a form on T*Q, then 

Lemma 6. Let X{t) be a curve in T*M with lf*^^{X{t)) = 'y{t). LetX{t) be 
an -horizontal lift of X{t) with U^'^{X{t)) = j{t). Then X{t) = A^(t) + 
X^{t), where is the curve in Ann(V) determined by 

(13) A^(t)(?y) = X{t){Tr,^ for any v G T^^t)Q, 

and X^ is a curve in Ann('H) satisfying 

T^(,) A^(t) = 0. 

Proof. First observe that since vr^ commutes with addition on respectively 
T*Q and T*M. Hence, both and are preserved under addition. In 
particular, the sum of two 7^^-horizontal curves are again H?. 

Let A^(t) be defined as in (jl3p . We will show that this is an ?^^-horizontal 
lift. Clearly ■K'^{X^{t)) = A(t), so aU we need to show is that X^{t) is Ti^- 
horizontal. By using the fact that we, at least for short time can write A« as 
t I— > 7r*a|;^(j) for some one-form a on Af, it is simple to verify that s^(A^, •) 
vanishes on V^. 

Using the conclusions from the previous two paragraphs, we obtain that 
)^{t) = X{t) - X^{t) is an 'H^-horizontal lift of 0^*^ \^{t). Let 5 be a section 
of Ann(?^), which, at least for short time, satisfy 5;|;y(t) = A(t). Let T be any 
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section of V, and let be a section of which projects to T. Inserting 
this into the symplectic form on T*Q, we have 

ci{>y,W) = dta{T)\^(t) - Warn)) - a(T^(,) T) = (t^(,) a)(T) = 0. 

Here we have used that 7 is ?^-horizontal. □ 

5.1.3. Proof of TheoremUl Let H be any smooth function on T*M and let 
H = H o tt'^. By definition, it is clear that iH is the Hamiltonian vector 
field of H. 

To prove (a), assume that X{t) is a curve in T*M along 7(t). Let 7 be 
a horizontal lift and let be defined as in (jl3p . By Lemma |4] (b), we 
know that £H and h^H coincides on Ann(V). The result then follows from 
Lemma [6l 

For the proof of (b), introduce the map <I>^ : T*Q — )• T*Q defined by 
<I>^ {p) = p pry . Obviously 

$^*(/i|u;-£pU') = 0, 

from LemmaU] (b), since pJZ{hqli^* , •) vanishes on vertical vectors. As 
such, if A is an integral curve of H with <I>^(A) = v{t) and '7r^(A) = A, then 
clearly WjV = Ohy Lemma [Gj where 7 is the projection of A to M. We 
again use Lemma |4](b), for the conclusion 

A = 7r\£H\r = H\x - vU A7^(7, •) = H\x - vIa ^;7^(7, •) 



Finally, we know that 7 = n^*^(A) is an ?^-horizontal hft of 7 = n'^*^^(A) 
from Lemma m (a). 

5.1.4. Representation in local coordinates. Recall that M is n dimensional, 
while Q has dimension n -\- v. Indices running from 1 to n will be in latin 
letters i, j, k, while indices running from 1 to i/, will be in greek letters k, A, fi. 

Let (x, U) be a local coordinate system on M, and let [/ be a neighborhood 
of 7r~^(C/) such that V trivialize over this neighborhood. Let T]^, . . . , T^ 
be vector fields on U forming a basis for V. We give T*U coordinates 
p = Yll=iPidxi and T*U coordinates p = Yli=i + ^^^-^ ?*kT*. Here, 

T* is defined by formula 

T:(/ia,j = o, t:(t^) = <5,^. 

Relative to these coordinates, 

n V n 

i=l K=l i=l 

The bundle is consequently spanned by and Write 
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Then ? = X]j=i ^ dpi, while 

n oo 

?' = Tr*dxi A (ioj + T* A dbi 

i=l K=l 

u I n n u V 

+ E ^« E ^ + E E A t;^ + ci^Ti A 

K=l \ i,j=l i=l /i=l -^i/i=l 

It is then simple to verify that T-L^ is spanned by 

V 

h^d,^ = hd,^ + ^ 6^rf,9fe^ and h'dp^ = da,, 

or 

V n V 

id., = hd^-Y, E ^-^^^1 - + E and = da, . 

K=l j = l K,fl=l 

5.2. Proof of Theorem [3]^b). We will use the notation of Section 15.11 
First, consider any curve t i-^ a{t) in Ann(D) which has a projection 7 in 
M. For sufficiently short time, define 5 as the ^{^-horizontal lift. Then 
by Lemma El this is in Ann(D). Hence, if w € rAnn(D) C TTM, then 
hjw G rAnn(I)). 

Let 1 1-^ A(t) be a characteristic of D with 7r^(A(t)) = A(t). We also write 
n'^*Q(A) = 7 and n'^*^(A) = 7. Finally, let v{t) = X{t) pry. We then have 
the following identities that X{t) must satisfy. 

For any Sq € Ann(Z));y((), we need to have 

?(A(t),vl3;^^^ "0) = 00(7) = 0, 

which is the horizontally requirement. Furthermore, let X, be any vector 
field on M. Then 

?(A, hW*^\x) = c^(A, o^**^x) - m n{x, 7) 

= c^(A, 0^**^X) - 7^(X, 7) = 
Finally, for any T € r(V), we have 

?'(A,0^*^,T) = -d^^(x,0'^'Q,T) = -j^v{t){T) =0 

5.3. Equations for a pendulum with varying length. Consider a pen- 
dulum in with mass m, under the infiuence of a force mA(cos (/>o, sin c/jq). 
The pendulum is connected to a massless rod or string whose length varies 
according to the function p{t). We give the pendulum coordinates 

p{t) (cos 6 {t),sm{t)). Newtons equation for this problem is then 

m{p - pO'^) {cos 9, sin e) + m{p9 + 2p^)(- sin 6*, cos 6*) 
= mAcos{9 — (/)o)(cos 9, sin 6*) — mAsm{9 — (j)o){— sin 9, cos 9). 
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The radial part of this equation is just the equation describing the constraint, 
while the other part is the essentially equation (llOp . 
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